Non-Markovian effects at the Fermi-edge singularity in quantum dots 
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We study electronic transport through a quantum dot in the Fermi-edge singularity regime, placing 
emphasis on its non-Markovian attributes. These are quantified by the behavior of current noise 
as well as trace-distance-based measure of non-Markovianity and found to be pronounced at low 
temperatures where the interplay of many-electron correlations and quantum coherence present in 
the system leads to significant quantum memory effects and non-Markovian dynamics. 

PACS numbers: 72.70.+m, 73.63.-b 



I. INTRODUCTION 

The Fermi-edge singularity (FES) is a phenomenon 
originating from the interaction of conduction electrons 
with localized perturbations and is characterized by a 
power-law divergence. It was first predicted theoretically 
for X-ray absorption in metals [l[ and later verified ex- 
perimentally @, H| . The theory then developed [l|, 0-0| 
has been used to describe other situations involving a 
similar Hamiltonian and leading to the same power-law 
divergence, such as resonant tunneling through localized 
levels 0. 

The FES in transport through quantum dots (QDs) oc- 
curs in the regime where the energy of the localised QD 
level (s) is comparable to the Fermi energy of one of the 
leads The change of occupation in the local level (s) 
during the tunneling process leads to sudden changes in 
the scattering potential and hence, to singular behavior of 
the current through the QD at resonance and power-law 
dependence of the current away from resonance [cl-fl2l]. 
As shown in recent experiments [ill, EH , current noise at 
such a FES displays characteristic behavior, including a 
deep minimum and a slight super-Poissonian peak, which 
cannot be accounted for by the Markovian theory. Inter- 
play of many-body correlations and quantum coherence 
is responsible for the large non-Markovian corrections in 
the FES singularity region [13| . 

This paper has two objectives: first, we wish to elabo- 
rate on the theoretical method used for the explanation 
of the experiment in Ref. [13[ and, furthermore, we ex- 
plicitly demonstrate the presence of quantum memory 
by the recently developed measure of the quantum non- 
Markovianity based on the trace-distance [Til [To| . It is 
organized as follows. In the next section the system under 
study is introduced, the Hamiltonian governing its evo- 
lution is described, including the scattering term which 
is responsible for the FES behavior, and the equation of 
motion for the reduced density matrix is formulated. In 



Sec. IIII1 the method for calculating the tunneling rates 
is outlined which allows the non-Markovian memory ker- 
nel to be found. Sec. IIVI contains the formulas for the 
mean current and non-Markovian noise, as well as the 
results. These include a comparison between Markovian 
and non-Markovian effects, a description of parameter- 
dependences of both current and noise, and remarks on 
fitting experimental data. Finally, next Sec. IVl explicitly 
shows the non-Markovianity of the time evolution by us- 
ing the trace-distance-based measure of memory. Sec. IVII 
then concludes the paper. 



II. THE SYSTEM AND ITS EQUATION OF 
MOTION 



The system under study inspired by the experiments 
[ill EH consists of a quantum dot with single spin- 
less level embedded between two leads of non-interacting 
spin-less electrons, meaning that there is a single channel 
of transport responsible for the current between the emit- 
ter and collector leads. This may be easily extended to 
the spin- full case and/or a multi-level QD and resulting 
multi-channel transport. The FES occurs on the emitter 
lead (henceforth referred to as the left lead) with much 
smaller tunneling rate, even at the FES threshold volt- 
age, than the collector (or right) lead, which means that 
tunneling from the emitter is much less probable than 
tunneling out of the QD into the collector. The chemical 
potential of the collector is well below the QD level, while 
the emitter's Fermi energy can be tuned in the resonance 
with the level. For a schematic of the setup, we refer the 
reader to Fig. 1 of Ref. [l3[ . 

The Hamiltonian of the system is H = Hqd + Hlq + 
Hrq + Hlt + Hut + Hx, where Hqd is the QD Hamil- 
tonian 



Hqd = eocftd, 



(1) 
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with eo denoting the energy of the single level in the QD 
and d' , d being the creation and annihilation operators 
of an electron in the dot. H a o with a — L,R are the 
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Hamiltonians of the left/right leads respectively, 

H a o = e k a c k a C k a > 



(2) 



where the creation and annihilation operators for elec- 
trons in the left/right leads are denoted by cj. , c ka and 
the corresponding energies are e ka ■ Tunneling from the 
left lead into the QD and to the right lead is governed by 
the Hamiltonians H a x (a = £,i?), 



(3) 



Lastly, the Fermi edge singularity occurs due to the inter- 
action of the electrons in the left lead with the electron 
in the QD; this interaction is described by the scattering 
Hamiltonian Hx, 



where the first term on the right-hand side describes the 
effect of the traced-out right lead [Iff, while the second 
term accounts for the evolution resulting from the rest 
of the Hamiltonian, i.e. Hl = Hqd + Hlo + -Hlt + Hx- 
The joint QD and left lead density matrix is written in 
the vector form in terms of the QD states (|0) and |1) 
denote the empty and occupied QD, respectively) a = 

(fOQ) "ll) fQlj °"lo) 



Following closely the procedure of Ref. [18|, the next 
step is separating Eq. §5§ into four equations for the four 
elements of a. Tracing out the electron states of the left 
lead in the two equations for the evolution of the diagonal 
elements djj (j — 0,1), allows us to find the evolution 
equations for the QD occupations (po(t) for the empty 
dot and p\{t) for an electron in the dot), 



H x = 



E 

fcr, ,k' 



dU=V x dU 



(4) 



In this system, couplings of the left and right leads to 
the QD state have vastly different nature and are, there- 
fore, described by different methods. The right lead is 
not affected by the FES Hamiltonian, the lead is effec- 
tively empty (very low chemical potential) and under the 
experimentally well-justified assumption of constant tun- 
nel rate -y R = 2irJ2k R l*fcii| 2( K e ~ £ k R ) ( we use units in 
which K = 1 throughout this paper) the transport from 
the QD to the right lead may be treated exactly following 
Ref. 

This allows to write the equation of motion for the evo- 
lution of the density matrix of the QD and the left lead 
in a mixed form analogous to the description of charge 
transport through a Coulomb blockaded double quantum 
dot in a dissipative environment in Refs. [I?], Ell- The 
equation reads 



da(t) 
dt 






1R 











—JR 














-7fl/2 





o 








-7fl/2 



a(t)-i[H Ll a{t)}, 



dpijt) 

dt 
dpojt) 

dt 



dpojt) 
dt ' 

= 7flPi(*) - 2 E tk L Im [ Tri ( c k L voi(t))} 



(6) 



fe L 



(5) 



The remaining two equations (hermitian conjugates of 
each other) govern the evolution of <7 i (and <7 10 ) enter- 
ing the second of Eq. ^ and couple to the diagonal ele- 
ments ajj . To close the equations for the QD occupations 
Pj(t) we perform physically- motivated perturbative de- 
coupling (in relatively small t^ L ) of the diagonal elements 
of the density matrix into djj(t) = Pj(t) €5 Rj, where 
% = exp ( - [H L0 + jV x - ml E fc , 4 L c kL ]/k B T)/Zj 
are the grand-canonical density matrices of the left lead 
at equilibrium (described by the chemical potential 
and temperature T) when the dot is empty (J = 0) or 
occupied ( j = 1). Now we find 



Im [Tr L (c kL a i{t))} = Im 



+Vx)t c kL e- lH ^d m {Q) 



dt'p (t')Re 



dt'pi(t')Re 



)'Tr L (e^° 

(^_ l£0 ) (t _ t ')^^ Tr 



J e i(H L0+Vx )(t-t') c iH LO (t-t')x o 



(J) 



which, after inserting into Eq. ([5]) and substituting t' — > t — t', leads to the non-Markovian equation of motion for the 
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occupations 



d J^ = - fdt'Mt-t^Re 
at In 



- zeo > t G Q (t') 
t 



7RPx(t)+ / dt' Pl (t-t')2Re 



e -^-* e °J t 'Gi(t / ) 



= - dt'j L (t')p {t-t')+ 7RPl (t)+ dt' 7 £(t')pi(*-*') + T(t) 
Jo Jo 



containing the initial-condition-dcpcndent inhomogeneous term 

T(t) = -2^t kh Im e-^f-^^YL (V^+^c^e-^'^tO)) 

and the FES Green functions defined as 

G (t)= £ ^ L t ki Tr L [cl^e^+^c^e-^Ro 

Gi(t)= tkJ^TrL [e^+^c^e-^c^R, 

k L ,k' 



•T(t) 



(8) 



(9) 



(10) 



r 



III. EVALUATION OF THE MEMORY KERNEL 



The evaluation of the Green functio ns (TTOl) is the cen- 
tral topic of the FES theory @, H 0, EMiioi a review 
see Ref. @). The standard procedure is to treat the 
Green functions in terms of open line Lj(t) and closed 
loop Dj(t) factors, such that Gj{t) = Lj(t)Dj(t). Fol- 
lowing Ref. [![, these terms are found to be 



D /i(t) 
Lo(t) 

Li(t) 



k L 

E< 



(11) 



with S being the phase shift associated with the scatter- 
ing effects of the occupied QD level described by Vx [1] , 
£o the ultraviolet cutoff, and np(e) denoting the Fermi- 
Dirac distribution of the left lead. 

At zero temperature, the Fermi distribution reduces to 
np(e) = 0(^l — e), where Q{x) is the Heaviside step func- 
tion. This considerably simplifies the problem of find- 
ing the functions Lj(t) and leads to the known time- 
dependence [J |i[ of the zero-temperature FES Green 
functions 



G 0/ i(t) = ^ r(l - -) e-^iTitta 



(12) 



Here, a = — — is the FES critical exponent (governing 
the power-law dependence typical for the FES), 7 £ = 
2ir J2k L \tk L | 2 ^( e — e k L ) is the non-interacting tunnel rate 
to the left lead stemming from the Hamiltonian ([3]), and 
T(a;) denotes the Euler gamma function. 



Knowing the Green functions, we also know the non- 
Markovian rates 7 [ b) (£) = 2Re e'^f^'^G^t) 

Their Laplace transforms needed for the evaluation of 
the stationary current and noise read 

7L (z;A) = 2ir(i--)r(a)^ 



x Im 



• + 3a(l-tA) 



(13) 



7 £(z;A)= 7L (z;-A), 



being parametrized by the dimensionless energy-distance 
from the FES edge A = 

The finite temperature result is obtained (as shown in 
Refs. 0, H H3l) by substituting time t in Eq. ^ by 



the temperature-dependent expression 



sinh(7rfc B Tt) 



The 



Tvk B T 

tunneling rate in the Laplace space jl(z; A) is then equal 
to 



7T V 7T I 



2nk B T 



x B 



1 



^(1-*A) 



27rk B T 



-,a 



(14) 



where B(x,y) is the Beta function. The second of Eq. 
(Q~3|) for the back-flow tunneling rate ^(z; A) carries over 
to the finite-temperature case. 

The non-Markovian generalized master equation 
(GME) for QD occupations p = (po,pi) T stemming from 
Eqs. ([8]) and (J6j> can now be rewritten in the matrix form 



dp{t) 
dt 



= f dt'W(t')p(t-t') + T(t), 
Jo 



(15) 
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with the memory kernel 



W(t) = 



~lL{t) 

7i(*) 



-7t(*)-7fl 



(16) 



and the inhomogeneity T(t) = (T(i), — T(i)) T . The ex- 
plicit form of this initial-condition-dependent term de- 
scribing the effect of initial correlations in the system 
is not necessary as the stationary properties (mean cur- 
rent and noise) do not depend on it, since it vanishes 
fast enough for large times. It also drops out from the 
expression for the trace distance, see below in Sec. [Vj 

For evaluation of the current and noise we need to ex- 
tend the GME by including the counting field(s) and ex- 
press it in the Laplace space [l?], [H| as 

zp({x}, z)-p({ X }, t = 0) = W({ X }, z)p({ X }, z)+T({ X }, z) 

with the modified kernel (for details of introduction of 
the counting fields see Refs. [HI, (HI) 



(18) 

where {%} denotes the pair of counting fields xl, XR a * 
the left/right junction respectively. 



IV. CURRENT AND NON-MARKOVIAN NOISE 

Following the scheme for calculating zero-frequency 
cumulants when given the Laplace transformed mem- 
ory kernel of a generalized master equation described in 
Refs. [13)13) the mean current and non-Markovian noise 
in the FES regime can be found. The kernel is given by 
Eq. (fT8"|) and the resulting stationary current reads 



1(A) = e- 



7ii7i(0; A) 



' 7Jt + 7 L (0;A)+ 7 £(0;A)- 
The non-Markovian Fano factor is given by 
2 7im (0;A) 



(19) 



F(A) = 1 - 
+ 2 7 ,r 



( 7 fl + 7L(0;A)+7i(0;A)) 2 

(77? + 7l (0; AM (0; A ) - 7l'(0; A) 7 l (0; A) 



(77? + 7l(0;A)+ 7 £(0;A)) 2 



where the prime denotes the z-derivative at 



7^(0; A) 



7i b) (*;A) 



(20) 



0, 



The Markovian noise 



z=0 



is given by the part of the expression without the deriva- 
tives, i.e. the first line of Eq. ([2D]) only. If, additionally, 
no back-flow is taken into account 7^(0; A) — > 0, E q. f2"0l 
reduces to the widely used Markovian expression jll| . 

In the following, the parameters a = 0.45 and 77? = 
3.4 ■ 10 10 s _1 measured in Ref. [ll[ are used, unless 
explicitly stated otherwise. The a-dependent prefactor 
7£r(l — ^)£o/7r can be used to control the height of 
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FIG. 1: a) Zero-temperature current as a function of A. b) 
Non-Markovian (black) and Markovian (red) Fano factor at 
zero temperature as a function of A. 



the current peak; in Figs. [TJ @] EJ and [5] it is taken such 
that 7l(0,0)/7t? = 0.09 at zero temperature. 

In Fig.[T]the current and the Fano factor near the Fermi 
edge singularity at zero-temperature are depicted as a 
function of the distance in energy from the singularity 
which is represented by the parameter A. Hence, the 
energy distance is given in units of half of the tunneling 
rate from the QD into the right (collector) lead. This 
dimensionless parameter is a more convenient theoretical 
measure of the energy difference than the applied voltage 
Vs d , since the relation of the actual change of the energy 
difference //£,— eo to the change in applied voltage involves 
a leverage factor [H,[ll| 77 = — £o)/e(VsD — Vq) (where 
Vo is the FES threshold voltage) which is often hard to 
determine experimentally. 

The current shows a typical FES behavior (Fig. Q]})), 
i.e. sharp growth on the low-energy side of the singu- 
larity and power-law decay on the high-energy side. In 
Fig. Hb), the Fano factor is plotted for both the Marko- 
vian approximation (red line) and the full non-Markovian 
result (black line) . Comparison of the two curves shows 
that near the singularity non-Markovian features cannot 
be neglected, since they change the qualitative and quan- 
titative noise features. Firstly, there is a significant drop 
of the Fano factor minimum on the high-energy side of 
the singularity. Secondly, on the low-energy side a new, 
super-Poissonian peak appears. Both these features can 
be regarded as signatures of strongly non-Markovian dy- 
namics. 

Fig. [2] depicts the dependence on the critical exponent. 
The black curves have been transferred from Fig. [TJ the 
blue ones correspond to a — 0.35, the red to a = 0.55 and 
the orange to a = 0.65. The parameters 7° have been 
adjusted for the three new curves so that the current 
peaks are of the same height as for the black lines. As 
expected, energy dependence of the current on the right 
side of the singularity changes according to / oc A~". 
As for the Fano factor, for higher critical exponent val- 
ues the minima are deeper and the curve on its high- 
energy side is steeper. Varying the critical exponent in 
the studied range has little effect on the Fano factor on 
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FIG. 2: Normalized current (a) and Fano factor (b) at zero 
temperature as a function of A for different values of a = 0.45 
(black), 0.35 (blue), 0.55 (red), and 0.65 (orange). 



the low-energy side of the singularity. 

The dependence on 7^ is shown in Fig. |3] The black 
curves have been again transferred from Fig. [TJ while 7^ 
is halved for the blue curve and doubled for the red curve. 
The variable A depends on 7^ itself, so the actual energy 
changes twice as fast for the black curve and four times 
as fast for the blue curve with respect to the red curve. 
As in Fig. [2] the parameters 7° have been adjusted in 
both cases so that the current peaks are of the same 
height. For the current, changing the right tunneling 
rate results primarily in stretching or squeezing of the 
peak (in energy), while the Fano factor undergoes large 
variations. Firstly, the minima and maxima are highly 7^ 
dependent; the dependence is more pronounced than the 
a dependence. Secondly, the speed at which the Fano 
factor approaches unity on the high-energy side of the 
minima is significantly varied. Contrarily to Fig. [21 the 
rise on the right side of the minimum is steeper when 
the minimum is shallower. Both of these features are 
important when fitting experimental data and may serve 
to check the correctness of the measured right tunneling 
rate tr. 

The temperature dependence is depicted in Fig. |4] 
The current peak is damped and broadened for grow- 
ing temperatures, as well as blue-shifted, in correspon- 
dence with experimental data and theoretical predictions. 
The features of the Fano factor are also damped, and the 
super-Poissonian peak disappears at high temperatures. 
In fact, it is the non-Markovian features becoming less 
pronounced; they lose any relevance at higher tempera- 
tures. This can be seen in Fig. [5] where the Markovian 
(red line) and non-Markovian (black line) Fano factors 
are compared for the same temperatures as in Fig. [4] 
The crossover temperature for the irrelevance of non- 
Markovian corrections depends on system parameters a 
and jr, roughly being of the order of 7^, see Sec. IVl 

Finally, let us briefly comment on some consequences 
of our results for the comparison with experiments. One 
of the methods used for finding the critical exponent a 
relies on collapsing of the current curves corresponding to 
different temperatures to a common curve daEl. This 
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FIG. 3: Normalized current (a) and Fano factor (b) at zero 
temperature as a function of A for different values of ir = 
3.4 • 10 10 s" 1 (black), 1.7 • 10 10 s" 1 (blue), and 6.8 ■ 10 10 s" 1 
(red). Note, that the parameter A stands for twice the actual 
energy for the red curve and half for the blue curve with 
respect to that of the black curve. 




FIG. 4: Current (a) and Fano factor (b) as a function of A 
at different temperatures T — K (black) line, 0.1 K (red), 
0.5 K (orange), and 1.0 K (green). 



is done by assuming the proportionality 

1 - a ir](V S D - V )e 



I oc Im 



27rfc B T 



B 



2i:k B T 



■ a 



( 21 ) 

where Vs d is the source-drain voltage responsible for the 
QD level energy shift, Vb is the FES threshold voltage 
corresponding to the resonance of the QD level with the 
left lead's chemical potential /iz, and ij is the leverage 
factor, such that r)(V SD - V )e = A^f. If Eq. dHJ) cor- 
rectly describes the measured I —V curves then rescaling 
the current / — > It = IT a and the source-drain voltage 
Vsd —> Vt — {Vsd — Voje/ksT which lifts any tempera- 
ture dependence from it, 



It oc Im 



B 



1 



ir/Vr 



-,a 



(22) 



will reduce all temperature-dependent curves to one for 
the correct critical exponent. 

Eq. (|2"Tj) is a twofold approximation. Firstly, the cur- 
rent, Eq. (|19[) . is a non-trivial function of the tunneling 
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FIG. 5: Markovian (red) and non-Markovian (black) Fano 
factor as a function of A at 0.1 K (a), 0.5 K (b) and 1.0 K 
(c). 




FIG. 6: Current in the FES regime (black) at zero tempera- 
ture compared to the approximate current (|21|l used in scaling 
at zero effective temperature (blue) and 0.19 K (red). 



rates, for which the proportionality should hold. It re- 
duces to I ~ e7i(0; A) for jr ^> 7l(0; A); the theoretical 
analysis presented in this paper requires this condition to 
be fulfilled, nonetheless, small deviations of the current 
from the left tunneling rate should be expected near the 
singularity. 

Secondly and more importantly, comparing Eq. (|21j) 
with Eq. (|14[) shows an additional real, temperature de- 
pendent factor in the Beta function 7fl/(47rfcsT). This 
factor stems from the first term on the right hand side of 
Eq. which is responsible for the effect of the right lead 
on the QD. It describes the damping of the coherence of 
the QD state and is typically neglected. Although this 
factor is negligible away from the singularity, it plays a 
significant role in reducing the FES peak. The discrep- 
ancy is usually handled by introducing an effective tem- 
perature ksT* = (fcsT) 2 + r| , where the temperature 
enhancement is wrongly attributed to a non-zero intrinsic 
linewidth of the QD state Ti ascribed to the hybridiza- 
tion of the QD electron state with the emitter (left) lead 
states. 

The current at zero temperature as described by 
Eq. (|19p is compared with the corresponding curve given 
by Eq. ([21]) at K and 0.19 K in Fig. H The effective- 
temperature treatment reproduces the current curves 
given by Eqs. (|19l) and (fl~4")l at relatively high temper- 
atures; it is then possible to collapse them to a common 
curve and find the critical exponent a with reasonable 
precision (if the leverage factor 7/ is known). The tem- 
perature limit of its applicability depends on the critical 
exponent and the tunneling rate of the right lead 7ij and 
is of the order of 0.1 K. At very low temperatures, col- 
lapsing different temperature curves does not yield the 
correct critical exponent. 



V. DIRECT MEASURE OF 
NON-MARKOVIANITY OF THE TIME 
EVOLUTION 

In this final section we address the non-Markovianity 
of the evolution directly in the time domain using the re- 
cently developed [lij measure of quantum memory based 
on the trace distance of two initial conditions, which 
was successfully experimentally tested in a designed 
quantum-optical experiment [151 ] . The measure is con- 
structed as follows: first, for any two density matrices px,2 
their trace distance is defined as D(p\, p 2 ) — \ Tr \p\— p 2 \ 
with 1 = V MA. In our case, we only deal with diago- 
nal elements of the density matrix p — diag(po,Pi), which 
together with the normalization condition po + p\ = 1 
yields for two probability distributions pi i2 

D( Pl ,p 2 ) = iflpg -pl\ + \p\-pj\) = \pl -pl\ = |E|. 

Now, the measure of non-Markovianity Af of a dynamical 
map inducing time evolution (here basically the evolution 
kernel W(t), Eq. pH) is defined as 0, GJ] Af(W) = 
max pi 2 ( ) Jq >0 dtZ)(pi(t), p 2 (i)), i.e. sum of stretches of 
the distance, where the distance is increasing, maximized 
over the pairs of initial conditions. The increase in the 
distance is a sign of back-flow of information into the 
system and is a hallmark of non-Markovian behavior [lij . 

We can find easily the time evolution of the quantity 
£(i) = Po(t) — Po(t) from Eq. (fTB"j) leading to the expres- 
sion in the Laplace space 

= s(t = o) = nt = 0) 

z - Wn(z) + W 12 {z) z + 7 fl + 7i(z) +7^(z)' 

(24) 

Obviously, E(i) (given as the backward Laplace trans- 
form of the above expression) is simply proportional to 
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FIG. 7: Memory measure E(t) (see the main text for details) 
as a function of the dimensionless time jut for various values 
of temperature k B T/-y R = (blue), 1/10 (purple; T = 0.03 
K), 1/3 (red; T = 0.09 K), and 1 (green; T = 0.26 K) at 
A = 0. 

the initial condition £(i = 0) — maximization over the 
initial condition in the measure AT(W) requires taking 
its maximal value £(i = 0) = 1. The measure uses in its 
definition the absolute value of £(t), which indeed may 
become negative, in order to insure the positivity of the 
distance. Straightforward inspection of the behavior of 
the distance reveals that the non-Markovianity measure 
can be obtained in our case by the same prescription ap- 
plied directly to £ and, in practice, this means that the 
non-Markovianity measure is basically equal to the depth 
of the E(t)-curve minimum, see Figs. [71 [HI 

In Fig. [7J we show the time evolution of £(<) at A = 
for varying values of temperature, while Fig. [8] depicts 
the A-dependence at zero temperature. Most impor- 
tantly, we see that the memory measure is generically 
non-zero (i.e., there is a finite dip below zero in the 
curves whose depth gives essentially the memory mea- 
sure) thus proving by a different method complementary 
to the noise the existence of the quantum memory in the 
system. In Fig. [Jj we observe an expected monotonous 
decrease of the non-Markovianity measure with increas- 
ing temperature fully in line with previous findings on 
noise shown in Fig. [SJ In Fig. [5] there is, however, a 
surprising non-monotonic dependence on A, namely, 
the curve for A = 1 does not exhibit any minimum and 
the measure assumes zero value. Moreover, Eq. (|2"4")l 
contains only a symmetric-in-A combination of rates 
and the S(t)-curves as well as the memory measure 
are therefore equal for ±A in disagreement with the 
obvious A-asymmetry of the Fano factor curves in 
Fig. [TJj). This may be caused by the neglect of the 



FIG. 8: Memory measure E(t) (see the main text for details) 
as a function of the dimensionless time ynt for various values 
of A = (blue), 1 (purple), 2 (red), and 3 (green) at zero 
temperature. 



off-diagonal matrix elements in our theory for short 
times (they are irrelevant in the stationary state and, 
thus, do not influence the noise). Their proper inclusion 
(though not obvious how to accomplish) might increase 
the measure by enlarging the space over which the 
maximization takes place, and potentially also break 
the symmetry in A which is now present. Despite of 
this issue, the global picture in Fig. [5J i.e. studying not 
just a single curve but a whole set of them, clearly re- 
veals the presence of the quantum memory in our system. 



VI. CONCLUSION 

We have studied electronic current transport through 
a quantum dot in the Fermi edge singularity regime and 
have shown that the associated current noise displays 
pronounced non-Markovian features at low temperatures. 
This is due to the interplay of many-electron correla- 
tions and quantum coherence present in the system which 
leads to the significance of quantum memory effects and 
strongly non-Markovian dynamics. The features include 
a pronounced deepening of the minimum in the Fano 
factor on the high-energy side of the singularity and an 
appearance of a super-Poissonian maximum on the low- 
energy side. 

The study of the current and noise curves parameter 
dependence shows that although changing the critical ex- 
ponent and the tunneling rate of the collector (right) lead 
have similar effects on the current curve, it leads to qual- 
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itatively different variations of the Fano factors energy 
dependence. Hence, noise measurements may be used to 
check the correctness of the experimentally found values 
of the critical exponent governing the FES power-law di- 
vergence. This is important since there is an ambiguity 
when fitting the current curves alone (due to problems 
in establishing the correct leverage factor between en- 
ergy and applied voltage) and we have also shown that 
the method of collapsing different-temperature curves to 
find the critical exponent is approximate and fails at very 
low temperatures. 

Finally, we have checked our noise findings by ap- 
plying a newly developed measure of quantum non- 
Markovianity to our problem. The results indeed con- 
firm the presence of quantum memory in the dynamics of 
the FES transport, although the application of the mea- 
sure contains slight inconsistencies pointing towards in- 



adequacy of theoretical description for short times. This 
could be a topic of further research alongside an obvious 
task of deriving the governing master equation from a 
systematic microscopic theory, which is an open nontriv- 
ial issue at the moment. 
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